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We consider the recent experiments on He^ bi-layers^, showing evidence for a quantum critical 
point (QCP) at which the first layer localizes. Using the Anderson lattice in two dimensions with 
the addition of a small dispersion of the f-fermion, we modelize the system of adsorbed He"^ layers. 
The first layer represents the f-fermions at the brink of localization while the second layer behaves 
as a free Fermi sea. We study the quantum critical regime of this system, evaluate the effective 
mass in the Fermi liquid phase and the coherence temperature and give a fit of the experiments 
and interpret its main features. Our model can serve as well as a predictive tool used for better 
determination of the experimental parameters. 



I. INTRODUCTION 

In the last fifteen years, an increasing body of ex- 
periemental results has revealed remarkable properties 
in heavy fermions, close to a zero temperature phase 
transition^ ^ ^. The standard laws governing the behavior 
of metallic conductors at very low temperature appeared 
to be violated in heavy fermions. Proximity to a QCP 
was early invoked, to explain the experiments^ ^ but so 
far, this wide body of observations remains a mystery 
and a challenging open problem. 

Recently, a new experimental set-up was explored, 
showing signs of quantum criticality of the same nature 
as for heavy fermions^, but in a rather different system. 
It consists of two layers of He^ fermions adsorbed on two 
layers of He^; those themselves adsorbed on a graphite 
substrate. Tlie histo ry of He^ films adsorbed on graphite 
is quite rictP^^^^^. A first layer of He^ has been ad- 
sorbed on graphite in two typical situations; on top of a 
compressed He^ solid of density 11.2nm~^ and on top of 
a deuterium layer of density 9.1nm~^. In both cases, a 
solidification of the top He^ layer is observed at a ratio 
of densities Ni/Nguh = 4/7. This "magic" number corre- 
sponds to a half filled super-lattice of unit cell ^/l x ^/7 
(see Fig. [T]), formed on top of the triangular substrate lat- 
tice. Specific heat measurements show that the effective 
mass increases by a factor of ten in the approach of the 
transition. The magnetic structure of the localized phase 
has been extensively studied. It is believed to be a spin 
liquid induced by ring-exchange^. The precise determi- 
nation of this spin liquid phase, and particularly whether 
it is massless or massive, and whether it has some ferro- 
magnetic component is still under debat^. Then a sec- 
ond and third He^ layers were adsorbed. The originality 
of the experiment^ is that it is the first time that, when 
the second layer arrives at promotion, the first layer has 
not yet solidified. Hence there is a regime in coverage 
where the two first layers hybridize while layer one sits 
on the brink of localization. 




substrate 

Experimental details can be found in Ref.^. We give 
here a rapid summary of the main findings of this work. 
The second layer arrives at promotion at a total coverage 
oi N = 6.3nm~^. From 6.3 to 9.2nm~^, a carachteris- 
tic temperature Tq is extracted, from the specific heat 
measurements, below which the fluid bilayer has Fermi 
liquid properties with an enhanced quasiparticle mass. 
Above To, a Curie law is observed, as if the first layer 
de-confines from the heavy Fermi liquid and behaves as 
a localized spin while the second one behaves as a Fermi 
liquid. It is quite difficult however, to separate quantita- 
tively the contribution of each layer in the heavy Fermi 
liquid phase. This characteristic temperature seems to 
vanish at a coverage Ncrit = 9.95nm~^, the so-called 
"critical coverage" by the experimentalists, with a power 
law 

Tcoh - S^-^ , (1) 

where S = \Ncrit - N\/ Ncrit • 

The effective mass is shown to increase by a factor of 
18 at = 9.0nm~^ and seems to diverge at Ncrit with 
a power law 

m/m* r^5 . (2) 
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Beyond Nc^ the first layer is fully localized at all tem- 
peratures investigated. However, NMR studies show that 
at Nj = 9.2nm~^ the magnetization starts to grow in a 
rather abrupt manner. It is not excluded that a first or- 
der ferromagnetic transition occurs for N > Nj but an 
experimental evidence for it is still not conclusive. The 
localized phase is believed to be a spin liquid ; a small 
"bump" in the heat capacity marks the onset of the spin 
liquid parameter. Experimentally it is evaluated to be 
of the order of J ~ 7 mK. Last, an activation gap is 
extracted from the heat capacity measurements. It de- 
creases with increasing coverages and there are indica- 
tions that it vanishes at a coverage lower than Ncrit- 



In this paper we give the details of the calcula- 
tions whose results have been announced in a previous 
Lettei^^. We apply the theory of the Kondo breakdown, 
previou sly intro duced for the study of QCP in heavy 
fermionJ^^n^E3 to the system of He^ bi-layers. The for- 
malism is identical to the one developed in^^. We use the 
Anderson lattice model with the addition of a dispersion 
of the f-fermions, to describe the system of He^ bi-layers. 
The first layer, in the brink of localization, forms the 
lattice of f-fermions. When the first layer localizes, the 
lattice is half-filled by construction. Strong hard core re- 
pulsion is taken into account by a short range Coulomb 
repulsion [/, with U ^ 2{)K^ in agreement with the early 
studies of bulk He^^. The top layer is modeled as a free 
Fermi gas. Hybridization between the two layers consists 
of hopping processes from layer one to layers two and 
vice verssli^. 



The paper is organized as follows. In section [H] we 
present the Anderson lattice model and derive the slave- 



boson effective Lagrangian. Section [Hi] is devoted to the 
evaluation of the bare parameters' dependence in cover- 
age. This is necessary if we want to confront our theory 



to the experimental data. We present in section [IV| the 
mean-field approximation. We show the presence of a 
QCP at T = corresponding to the Mott localisation 
of He^ first layer's fermions. In particular, a peculiar 
behavior of the effective hybridization explains the ap- 
parent occurence of two QCPs in the experimental data. 
We then study the fluctuations in section[V|discussing the 
critical regime and computing the effective mass and the 
coherence temperature in an intermediate energy regime 
corresponding to a dynamical exponent z = 3. We con- 
clude in section [VIl| with our main result and give a crit- 
icism of our work. Some technical details are presented 
in the appendices. Appendix A shows the calculation of 
the integrals used at the mean-field approximation. In 
Appendix B, we give the details of the evaluation of the 
fermionic contribution to the corrections of scaling of the 
holon mass and discuss the stability of the QCP. Finally, 
in appendix C, we derive an expression of the free energy 
starting from the Luttinger-Ward formula. 



II. THE MODEL 



Our starting point is the Anderson lattice model 



^ = E [4 + (^^0 - ha 



(3) 



Here (i, j) refers to nearest neighbour sites created by 
the Graphite's corrugate potential, a is a spin index, 
ficrificr) ^rc crcatiou (annihilation) operators for the first 
layer's fermions, cJ^(Qcr) are creation (annihilation) op- 
erators for the second layer's fermions. tij = t is the 
c-fermion's hopping, t^j = at is the f-fermion's hopping 
term, V is the hybridization between the two layers, Eq 
is the energy level of the f-fermions and /i is the chem- 



ical potential. n/,i = ^^fj^fia and = Y.^c\^Ci^ 
are the operators describing the particle number of each 
layer's fermions. U and Ui are respectively the intra- and 
inter-layer Coulomb repulsion. The model is studied in 
the limit of very large on-site repulsion U. We expect to 
have a coverage dependant hopping parameter, t = t{N) 
as well as hybridization V{N). Furthermore, we have 
/71 <C /7, but we keep the inter-layer interaction term for 
now. 

Super-exchange terms can be generated by a second or- 
der expansion in large U/{at) and Ui/{at). The Hamil- 
tonian is then written 



n 



where J = 2{atf/U,Ji = 2{atf /Ui and S/ = 
f3 fa^apfp ^P^^ Operator with a the Pauli ma- 

trix. RKKY interaction, mediated by the conduction 
electrons, as well as various ring exchange parameters, 
studied in^^ can be included in the J-term. 

One key approximation of this work is that we consider 
that at the edge of localization, the f-fermions are half- 
filled. This means that the f-fermion somehow form their 
"own" lattice as the coverage increases, so that when the 
localization occurs, we are at half filling. This approxi- 
mation is necessary if we want to attribute the observed 
increase of the effective mass to strong correlations com- 
ing from Mott physics. However, we don't have a micro- 
scopic justification for it; only the coherence of the find- 
ings of this approach can justify it. The on-site Coulomb 
repulsion U is very large ( ~ 20 K)^ leading to strong 
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correlations effects. In the limit ^ oo, there is a con- 
straint of no double occupancy which we account for us- 
ing Coleman's slave bosoiP^, decoupling the f-fermion's 
creation operator at each site "i" as 



(5) 



where p , the creation operator of the so-called "spinous" 
and 6^, the one of the holons, obey the local constraint 
Scr fiafio- bjbi = 1. Upon the transformation ([s]), the 
slave boson drops of all bi-linear products of fields at the 
same site. 

The constraint is taken into account in a Lagrangian 
formulation through a Lagrange multiplier A. The effec- 
tive lagrangian is then 

i ia 

+ ^ J (S/,,.S/,, - n,n,/4) + JiS/,,.Sc,, , (6) 



where Cf^i = Eq — ja -\- Xi is the renormalized f-band's 
chemical potential. 

The short range magnetic interaction and the 
induced Kondo interaction are decoupled using 
Hubbard- S tr at anovich transformations : JSf^i.Sfj 
^ijflfjcT-\(l)ij\^/J and JiS/,i.Sc,^- ^ crjlci^-\ai\^/Ji. 

The Lagrangian becomes now 



JO- 



E 4idr + m + E {fUVh + CTi)Ci. + h.c) 

i ia 

Y.(h + w.\V{j,)-Y.\<Pi,\Vj (7) 



We assume that condenses in a uniform spin liquid 
phase, i.e. = ^o- It renormalizes the dispersion of 

the spinon band and ensures the breakdown of the Kondo 
effect^. It is shown to stay r oug hly constant through the 
phase diagram = /3t ^ fi^. cfi merely renormalizes 
the effective hybridization Vhi. 



III. THE PARAMETERS 

Before going further, we need to evaluate the depen- 
dance of the bare parameters in coverage in order to fit 
the experimental data. 



The height of the layers is taken from the study by 
Roger et al}'^: the first ^He layer's height is ^ 2.02 A 
while the others' one is ~ 2.85 A (see Fig. ^. From the 
experiment^, the density of He^ layers is 9.2 nm~^ while 
the one of the first He^ layer is A/"! = 6.3 nm~'^ . The total 



^He bilayer 



12.851 




FIG. 2: ^He solid layers on top of the triangular lattice of the 
substrate; we show here the various heights of the layers one 
compared to the other. 



coverage is defined as 



Nf = 7Vi(l-n5) 



(8) 



where Nf and Nc are respectively the coverage ( in 
nm~'^) of the first and second layers and is the number 
of holons per site. At the transition, we have 

Nf/Ni = 1 (9) 

which accounts for the fact that at the transition, the 
f-fermions are in a 1/2 filled lattice. This means that the 
number of holons rib vanishes at the QCP. Away from the 
QCP, the number of holons is allowed to fiuctuate freely 
and its value is determined self-consistently. 

The parameter J is extracted from the experimentm 
J = 7mK. 



The evaluation of the bandwidth, D = 2t = 7r/m, for 
each layer of He^ is based on an analysis in Pricaupenko 
and Treiner^^ where the kinetic energy of liquid He^ con- 
tains a density dependant effective mass : 



2m* 



2m 



P_ 

Pc 



(10) 



where p = 3/{4:hc)n{A~'^) is the average density inside a 
sphere of radius he = 2. 63 A and pc = 0.04A~^. 
We have then 

Df = D {1- 0.07 Nff 
Dc = D {1 - 0.07 {N - Nf))^ , 

where D is the bandwidth of ^He in the bulk. 
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At half filling, Nf = A^i, the mean kinetic energy 
Ekinj equals the bandwidth Dj. We have 



kin,f 



2m 



}Ni Jo 



(27r)2 



16m 



(11) 



/ 



where kp = tt/v^TVi {y/Ni/2 is the average radius of a 
particle in the first layer). 
We find then at half-filling 

Ecj ^ 0.52 K 
Df = 0.30D 
Dc = 0.62D , 

thus, D ^ 1.73K, Df ^ 0.57 K and ^ 1.18 K 
which gives a value a = 0.54 for the ratio between the 
bandwidths. This value is relatively high compared with 
the typical values for rare earth compounds for which 
a = 0.1. 

A word has to be said at this stage : we have considered 
the spherical dispersion of the free fermions 



2m 



2m 



for which the density of states (DOS), defined by ^ 
p{e)de^ is constant 



27T 



However, as emphasized in the introduction, the first 
layer solidifies into a triangular lattice. For a triangu- 
lar lattice tight-banding band structure, the dispersion is 
given by 

Ck = —2t ^cos (kx) + 2 cos {kx/2) cos {V3ky/2)^ . 

The Fermi surface for fermions in a triangular lattice 
is no longer circular at each filling, but we can consider 
that these deviations are benin in the range of coverage 
studied in our case, in particular very close to the QCP. 

Fig. ([3| shows the Fermi surface of the f-fermions 
at two different coverages : (a) J = 0.10 for which 
Df ^ 0.b2K,ef ^ 0.13K and (b) S = 0.15 for which 
Df ^ O.QOK^Cf ^ 0.54i^r. In the latter case, the Fermi 
surface deviates around the circular Fermi surface for free 
fermions. 

The approximation of constant DOS can still hold and 
this can be seen indeed by considering the DOS profile 
for the triangular lattice case shown in Fig. |4] The 
hatched region marks the energy scales of our model, 
and we see that we are far from the Van Hove singularity, 
and we can approximate the DOS by a constant one. 





(a) 



(b) 



FIG. 3: The Fermi surface at unit coverage : (a) 5 = 0.10 
and (b) S = 0.15. We see that the Fermi surface in the former 
is still circular while it experiences, for the second one, small 
deviation from the circular case. 




FIG. 4: Density of states for a triangular lattice tight-banding 
structure. Characteristic energy scales in our model lie within 
the hatched region. 



The chemical potential /i is defined by the filling of the 
second layer 



6.3 



pode 



We get directly 



^ ^ 6.3 ^ 



(12) 



£^0 is identified as the difference between the potential 
energies of the two layers^^. Each layer experiences two 
kinds of interaction : 

• Van der Waals interaction with the grafoil substrate 

V^iz) = (4C|/27Z)2) 1/^9 _ Cs/z^ , (13) 

where D = 192 is the well depth of the po- 
tential and Cs = 2092 KA^ is the Van der Waals 
constant^^, and 
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• the Bernardes-Lennard Jones interaction between 
two He particles 



VLj{z) = Ae{{a/rr-{a/rf) 



(14) 



with e = 10.2 and a = 2.56 A is the hard core 
radiu^^. Thus, for the layer L^, the potential en- 
ergy writes 



with 



Ef^c = Vs{zi) ^Vi , 



Vi=TT^pj j rdrVLj(r) , 



(15) 



(16) 



where pj is the density of each layer and "j" is the 
layer's index. 

the chemical potential now reads 

Eo = Ef-E,. (17) 

We denote (see Fig. [2| ri, r2, rs and respectively 
the distances of the first, second, third and fourth layers 
to the graphite center. We have 



2.2A , Z2 = 5.03A , Z3 = 7.9A , z^ 



10.57 A . 

(18) 

Applying ([l3| we get Vsizs) = -4.21 K and Vsiz^) = 
— 1.77 K. These orders of magnitude are quite big com- 
pared to the typical scale of a few mK for this system. 
Its order of magnitude is in accordance witlP^. 

We turn now to the Lennard- Jones potential. We sum 
up (14) for all two body interaction in all the layers. We 



get for the first layer or "f"-fermions 



Vf 



Pi 



Ps 



rdrVLj{r) + 
rdrVLj{r) + p4 



P2 



rdrVLj{r) 



rdrVLj{r) 



,(19) 



with 



pi 0.092 , p2 = 0.092 



5.7, 



2.85 , 



PS = 10-^Ni{l - Ub) , P4 = 10-^[N - Ni{l - Ub) 

ri = 2.85 . 



3 -1/2 



For the second layer or "c" fermions, we get 

Pi / rdrVLj{r)^p2 / rdrVLj{r) 



Vc = TT 



Ps / rdrVLj{r) ^ p4: / rdrVLj{r) 



,(20) 



with 



pi = 0.092 , p2 = 0.092 , 



8.55 



5.7 



PS = 10-^Ni{l - Ub) , p4 = 10-^[N - 7Vi(l - Ub) 
rl = 2.85 , 



4 _ -1/2 
— P4 



The values of rl^ are now in A. 




dx dy 
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We finally get 

Vf = -3.87- 1.08A^-6.75n5 + /3 (21) 
with Is = -0.196(1 -n5)^[11.2- 0.3(1 -n^)^] . 

and 

Vc = -7.69 + 6.75n5 + J4 (22) 
with J4 = -7rlO"^(7V- 6.25(1 -715)^) 

[0.29 -3.210"^(7V- 6.25(1 -715)^)] . 

Eq now reads 

^0 = 1.65 - 1.0717V - lO-^TV^ - 13.5n6 + 
-2.25(1-715)^ 

+0.059(1 - Ub)^ (23) 



The last parameter, and the most crucial in fact is the 
hybridization V. It is defined as the hopping strength 
between the two layers. We can have an estimate of V 
using equation (ITo| to get the same dependence as in ( 11 ) 



V = Vo{l - 0.07 (TV - 7V/))(1 - 0.077V/). 

Here, Vq is proportional to the overlap between the 
ground state wave functions of the two layers, i.e. 
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ti2-SV I dz<!>i{z)<!>2{z) 
where is the ground state wave function of layer i and 

6V = Vs{z4)-Vs{z3). 

The latter is taken as a Slater determinant of single 
particle states which writes, assuming translational 
invariance parallel to the surfacd^^^ : 

^?(r) = ^^H^) exp [i{k^x + kyy)] 
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Density functional models show a Lorantzian-like profiles 
for the density of each layer along the z-direction^^ : 



Pi l</>i(^)P 
FronFl, we have : 



ForLi: a = 1.7 6 = 0.115 
For L2: a = 3.42 h = 0.32 

We find then Vq ~ 0.6i^ consistent with the value 

P 




FIG. 5: Sketch of layers' density profile. The hybridization is 
estimated from the overlap between the wave functions of the 
two layers. 

obtained in a previous stud^EU. 

As said before, the hybridization V is actually a crucial 
parameter. Indeed, the me an fie ld value for the QCP 
reads J/t — exp [£^o^c/^^] ^^^^^ we see then that any 
small variation in the dependance of V on coverage has 
an exponential impact on the position of the QCP. We 
consider V thus as a fitting parameter that will tune the 
position of the QCP. 

We have taken 

V = Fo(l - 0.07 {N - Nf)){l - 0.07 AT/) + + ¥28'^, 

where Vq^ Vi and V2 are adjusted to fit the experimental 
data and 5 = {Merit - N)/ Merit- We used Vq = 1.55 AT, 
Vi = 15.9 and V2 = -4.5 . 



IV. MEAN-FIELD THEORY 

At the mean field level, we make a uniform and static 
approximation for the holon field and the Lagrange mul- 
tiplier. The free energy writes then 

Fmf =-2TEfc,.,..^,± ln(-Gii(zc^n,k)) 

+A(62 - 1) (24) 

where cOn is the fermionic Matsubara frequency and 
G^^ = icjn — £^k±, with 



k± 2 



±A/(ek 



4V252 



In the obove, ek is the dispersion of the conduction 
electrons, = {aP + P)e\^ + ^/ is the spinon dispersion 
and E\^± the dispersion of the renormalized upper (+) 
and lower (-) bands (See Figj6| The former derives from 
the c-fermions with weak f character whereas the latter 
derives from the f-fermions with weak c character. 




Renormalized 
bands 



FIG. 6: Sketch of the different dispersions: the lower band 
is the dispersion of the spinons and the higher band is the 
dispersion of the conduction electrons. The Kondo gap A is 
defined as the difference between the chemical potential and 
the upper band. 



Minimizing (24) with respect to the holon field b and 



the Lagrange multiplier A, one gets the following mean- 
field equations 



T E Gfj + b^ = l, 



(25) 



where 



Gfc 



{icJri - £^k+)(^^n - £^k-) 

Vb 

{icJn - £^k+)(^^n - £^k-) 
VbPfc 



(26) 



These equations are solved in the case of a linearized 
dispersion bandwidth at zero temperature (T = 0). The 
summation over (k, cj^) is performed anatically and is 
given in Appendix A. The set of resulting equations is 
then solved numerically. 

FigjT] shows the plot of the order parameter, defined 
as the effective hybridization Vb^ and the "Kondo gap" 
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A, defined as the energy difference between tiie ciiemical 
potential and tiie upper band (See Fig. [6|, as a function 

0^8=l-N/N,r^t• 

In our model, the Kondo gap is identified with the 
activation energy observed experimentally in the specific 
heat. We have two bands in the model, one for the 
spinous and one for the conduction electrons. At very 
low hybridization, when the bands just start to hybridize, 
there is no energy difference between the lower and the 
upper band. As the hybridization grows, the upper band 
becomes empty and an activation gap opens. We see on 
Fig. [7|that the gap closes at the very vicinity of the QCP. 

The set of mean-field equations shows a QCP where 
6^0, the so-called Kondo Breakdown (KB) QCP, which 
implies that the spinous experience a Mott transition and 
their band is half-filled. We observe that Vh goes to 
zero, before the experimentally observed QCP occurs, at 
a unit coverage 5 ^ 0.063. This constitutes one main 
finding of this paper. The localization occurs before the 
experimental QCP is reached. Our interpretation is that 
first, the experimental QCP is evaluated by extrapolating 
to zero temperature the power laws for the effective mass 
and the coherence temperature; second, a key feature of 
the model is that the hybridization is strong, compared to 
the other parameters (it is of the order of the bandwidth) , 
hence the falling down of the order parameter close to the 
transition is very abrupt. 

This fact is illustrated in Fig. [7] where we see that the 
order parameter's behavior has two regimes: it starts to 
grow very quickely at the QCP then reaches, at the "el- 
bow" , a regime of strong hybridization. The behavior of 
the order parameter is governed by the relative strength 
of the bare hybridization V compared to the other en- 
ergies of the model. The former is already big at the 
QCP, Vc ~ 1.63i^, thus the slope of the effective hy- 
bridization is steep in the hybridized phase. The sharp 
change corresponds to the emptying of the upper band, 
the same point at which the opening of the Kondo gap 
occurs. This point is situated after the real QCP, in the 
hybridized phase, because when the localization occurs, 
the f-band is half-filled and the upper band is constrained 
to sit below the chemical potential and is thus occupied. 
The vanishing of the Kondo gap before the QCP is ob- 
served experimentally if we identify it as the activation 
gap extracted from the thermodynamic measurements of 
Neumann et al} 

We can make the same construction as the experimen- 
talists, by extrapolating the order parameter in the high 
energy regime to zero temperature. We find an addi- 
tional QCP that we identify with the "experimental" one. 
This gives an explanation of the mysterious presence of 
two QCPs in this system; the magnetization starts to 
grow at the physical QCP, before the experimental one 
is reached. Indeed, as soon as the first layer localizes, one 
expects the static magnetic susceptibility to grow quickly 
since the spin liquid parameter is small J ~ ImK. Note 
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Effective Hybridization Vb 
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0.05| 0.1 

fexp.QCP real QCP 6 



0.15 



0.2 



FIG. 7: Mean- field phase diagram for the Anderson lattice 
model in D = 2 applied to He^ bilayers^^. Following^ ^ = 1 — 
N/Ncrit with Ncrit — 9.95nm~^ . The effective hybridization 
(red online) drops suddenly at ^ ^ 0.063, indicating the real 
QCP. The experimental QCP is obtained by extrapolation of 
Vb to zero (E=0). The Kondo gap A (black online) vanishes 
before the real QCI^. 



that the distance in coverage between the two QCPs is 
in agreement with the experimental data. 




0.25 



FIG. 8: The number of bosons Ub (green online) and of the 
f-fermions (red online) per site in the system as well as their 
sum (black). The local constraint of no double occupancy is 
preserved throughout the explored coverage range. 

In the Figure |8] we have plotted directly the number of 
holons in the hybridized phase, as given from our mean- 
field theory. The number of holons determines the num- 
ber of holes in the first layer as compared to the value 
at half-filling. We can see that although the order of 
magnitude is correct close to the QCP, far away from it 
we obtain some values of too big from what is ob- 
served experimentally. In particular, it is believed that 
close to the coverage corresponding to the promotion of 
the second layer, the number of holons should decrease 
so that the number of f-fermions in the first layer should 
be again close to half-filling. We don't observe any hint 
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of this decreasing. It shows that the domain of vahdity 
of our model is close to the QCP. F ar away from it, we 
miss the physics of exhaustiorP^l^, where there are not 
enough free fermions in the second layer to Kondo screen 
the many f-fermions in the first layer. 



V. FLUCTUATIONS 

In what follows we will be interested in fitting the ex- 
perimental data. We identify the regime of critical fluc- 
tuations experimentally accessible with the higher energy 
regime of the order parameter (see Fig. [t]). Within 
our theory, we are situated in the intermediate regime 
around the Kondo breakdown QCP, i.e. the regime for 
which the dynamical exponent z = 3. We refer the 
reader to pr evious stu dies of the Kondo breakdown for 
more details^^EHIlIII'^ Xo give a small summary of the 
situation (see Fig. |9|, the main finding of the Kondo 
breakdown QCP is its multi-scale character. There ex- 
ists an energy scale differentiating two regimes. In 
the low temperature regime we have the dynamical expo- 
nent z = 2^^, with no damping. In the high temperature 
regime, we have the exponent z = 3 an the bosonic mode 
corresponding to the fluctuations of the order parameter 
is over-damped by the particle-hole continuum. In this 
paper we focus on the z = 3 regime, arguing that is 
very small in this system. Indeed, from the theory (see 



TO 



Regime II 



Hence we obtain 



= 8.10"^i^ 



K 



(27) 



which is a too small energy scale to be accessible experi- 
mentally for this set-up. 

The holon propagator in the intermediate regime (z=3) 
reads: 



D-ifeJ^„)=Do"' 



7l»n| 

a'q 



(28) 



with = 4kj,/{poV^), 7 = mV^Do/{7rvF) a' = 
b'^a + J/t, po = mc/(27r) is the c-fermions density of 
states and ^ is the correlation length, associated with 
the fluctuations of b, given by 
is the holon mass at T = 0. 



Dr. ^m. ^ where rub 



A. The Holon mass 

The static part of the holon mass is evaluated by dif- 
ferentiating twice the mean- field energy ( 24 ) with respect 
to the holon field b given the contraints (25). One finds 



rub 



26T^ 



dG 



ff 



k,a; 



db 



dP 



db 



(29) 



The summation over (k,a;) is evaluated anatically for a 
linearized dispersion bandwidth at T = and the result 
is given in Appendix A. 




QCP 



The temperature dependence of the holon mass is com- 
puted by evaluating the corrections to scaling to the bo- 
son propagator. There are two types of corrections to 
scaling. One contribution is the renormalization of the 
boson propagator coming to their coupling to the fermion 
loops. 



FIG. 9: Schematic for the Kondo breakdown QCP in the 
Anderson latticj^. On the left, where the holons are not 
condensed is the localized phase. On the right is the heavy 
Fermi phase. The QCP is multi-scale ; for T < E* , the 
dynamical exponent is z — 2 and for T > E* it is 2; = 3. 

for exampl^ we know that E* O.l^q* / qfTx, with 
g'* the mis-match of the two Fermi surfaces at the QCP. 
Here Tk can be taken as the typical energy scale of the 
system which is typically of the order of Tk = 100 mK. 
At the QCP, we evaluate q* /hp which is 

q* /kp = 1 - kc/kp , 

= 1 - (6.3/9.95)^/^ , 
= 0.2 . 




This type contribution was first evaluated close to a QCP 
in Ref.^. We first note that the two diagrams are propor- 
tional: li^a^ = a' n^^^ and that there is no corresponding 
vertex insertion at the first order. Hence, although the 
QCP occurs in the charge channel, we have no cancel- 
lation of this set of diagrams. This is in deep contrast 
to what occurs close to a ferromagnetic QCP or in the 
theory of non analytic corrections to the Landau Fermi 
liquid, where this set of diagrams cancels in the charge 
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channeP^. This type of diagram is known to be dan- 
gerous, and carries a minus sign, which destabihzes the 
fixed point. The diagram for U.^^ is computed in the in- 
termediate energy regime with the dynamical exponent 
z = 3. 

On the other hand, we have the direct mass renormal- 
ization coming from the standard (/)^-type corrections to 
scahng, which has to opposite effect of stabihzing the 
fixed point 



b b 



Here comes from the quartic term of the holon action 
derived from in a Ginzburg-Landau approach and 
contains the ferromagnetic short range correlations J ; 
we find ^4 = - J/4 + V^/{2a'^D^). We have as weh as 
the correction to the boson mass coming from the gauge 
fluctuation, which stabilizes as well the QCP, but is sub- 
dominant compared to the two previous ones. 



where Fm f is the free energy at the mean- field ( 24 ) and 



D{q^ VLn) is the full propagator of the holons. Note that 
we have neglected the role of the gauge fields in this for- 
mulation, because the re-normalization of the effective 
mass is to be evaluated inside the ordered phase where 
the gauge fields are gapped through the Higgs mecha- 
nism. At the mean-field level, the system consists of the 
upper and lower bands, we get then : 



(32) 



where is the density of states at the Fermi surface 

of the upper (lower) band given by : 



P± = Po 



k± 



|£;±=o 



The calculation is done in Appendix [C] 
The effective mass reads directly 

7^ 



m 



2^(p+ +/>_) + 



Summing the dominant diagrams (see Appendix [B| 
yields a logarithmic correction to scaling 



m5(T) = m5(T = 0) + C TLogT 



(30) 



where C had to be adjusted to C = 7.5 10 ^ to fit the 
data, while the analytic evaluation gives 



C 



(1 + aO 



DJ 



The balance of the two contributions in favor of the ^4 
coupling ensures the stabilty of the fixed point. 
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B. The effective mass 

The effective mass m* is determined from the free en- 
ergy of the system by : 

F = m 

6 

The free energy is evaluated using the Luttinger-Ward 
functional^^ 

F = Fmf+T/2 f (fq/{2Trflog [D-\q,nn)], 

n ^ 

(31) 



FIG. 10: The inverse effective mass m/m* in the Anderson 
lattice model for the He^ bi-layers^^. The dots are experi- 
mental data from^ . The fitting parameters for this model are 
detailed in the texPl. 



The result for m* is shown in Fig. 10 where it is com- 
pared to the results of experiment^. We see that the in- 
verse effective mass follows the same behaviour as the or- 
der parameter and vanishes at the theoritical QCP. Here 
again, if we extrapolate the high energy regime down to 
zero temperature, we can identify a fictions point where 
the effective mass could vanish, if it has not its peculiar 
behaviour into two regimes. This extrapolation is linear 
and follows closely the one found by the experimentalists. 
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C. The coherence temperature 

The coherence temperature is defined by the cross-over 
condition 

mb{Tcoh) = 0, 

where mi){T) is the temperature dependant holon mass, 
given in (30). 



The equation is solved numerically, using the re- 
sults found for the order parameter b in the previous 



section , and the result is plotted in Fig. 11 
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FIG. 11: The coherence temperature Tcoh in the Anderson 
lattice model for the He^ bi-layers. The dots are experimental 
data fronl^. The fitting parameters for this model are detailed 
in the text"*^^. 



The coherence temperature has the same qualitative 
behavior : it vanishes at the real QCP and we can 
extrapolate its high energy regime down to zero temper- 
ature closely to a quadratic power law in unit coverage 
8. 

In fact, the exponents of the effective mass and the 
coherence temperature can be understood in a simple 
way. For 2: = 3 theories in the Fermi liquid phase, the 
effective mass goes like the correlation length^ m/m* ~ 
From the dispersion of the boson mode we see that 
^ y^m5 ~ h. Now the coherence temperature goes 
like . In the regime where h varies linearly with the 
coverage n we thus get 



m/rri' ^ est 



Tcoh ^ {est - nf 



(33) 



VI. DISCUSSION 



One of the main general observation one gets from 
the experimental data is the asymmetry of the phase 
diagram, as far as the quantum fluctuations are con- 
cerned. Indeed the increase of the effective mass appears 



only from the right of the phase diagram which corre- 
sponds to low doping (see Fig. 10). From the left of 



the phase diagram the fluctuations seem to be frozen out. 

Another observation is the quasi-absence of quantum 
critical (QC) regime in temperature for this system, 
unlike for the heavy fermions. Indeed a Curie law for the 
spin susceptibility is observed at very low temperatures 
in the localized phase and directly above Tcoh in the 
hybridized phase, indicating that the system very quickly 
goes into a regime of free spins, hence missing the usual 
quantum critical regime typical of QCP. 

The key to understanding these two observations is 
that in this system the energy scales are completely 
different form the ones that appear in heavy fermion 
systems. 



50 mK 




Heavy Fermi 
Liquid 



V 



FIG. 12: Asymmetry of the phase diagram. Above T* the 
entropy RLn2 is released. There are two ways of quenching 
the entropy, first through the formation of the heavy Fermi 
liquid phase where the hybridization is non zero ( on the right 
of the phase diagram), second through the formation of the 
spin liquid (on the left of the phase diagram) . 



The Curie law is observed when the entropy Rhi2 is 
released, above a characteristic temperature T*. In our 
model, two mechanisms are responsible for quenching the 
entropy, namely the formation of the spin liquid and of 
the heavy Fermi liquid. T* is thus determined by the 
relative strength of these two mechanisms. Technically, 
T* is by the first irrelevant operator of the theory. We see 
in Figure [12] that on the left side of the phase diagram, the 
main quenching mechanism corresponds to the formation 
of the spin liquid, while on the right side of the phase 
diagram, the two mechanisms coincide and are roughly of 
the same strength. The asymmetry of the phase diagram 
can thus be accounted for, in this model, by the fact that 
on the localized side (left side) the spinous' bandwidth, 
which determines the scale of the the formation of the 



11 



spin liquid, is typically given by the value of the exchange 
parameter J 7 mK. Alternatively, in the hybridized 
phase, the bandwidth of the spinous is enlarged, due to 
the holon fluctuations Df = J -\- UbaD. This increase 
of the bandwidth in the hybridized phase is typical of a 
slave-boson description of a Mott transitiorP^. 




FIG. 13: The effective bandwidth of the spinons Df — J -\- 
aDuh and the experimental characteristic temperature Tq. 

In the hybridized phase the coincidence, within the 
experimental uncertainties (between 5 and 10 mK^), in 
energy between the cross-over coherence temperature 
and the effective mass of the spinons (See Fig. 



13) 



explains that the quantum critical regime is quenched, 
the free spin behavior being admittedly quickly observed 
above the temperature which delimits the upper-critical 
regime. 



VII. CONCLUSIONS 

In the present article, we give the details of calcula- 
tion whose results have been presented in a previous 
Letteil^. The system studied, He^ bilayers, is one of 
the simplest physical ones, with negligible spin-orbit 
interaction and no crystal-field interactions, to show 
quantum criticality (QC) similar to the one observed in 
complicated intermetallic heavy fermions compounds. 



Using the Kondo-BreakdowrP^l^^^ scenario of an 
itinerant QCP, we examine the possible origine of the QC 
observed experimentally as fluctuations of an effective 
hybridization. The theoritical model is an extended 
version of the Anderson lattice model with a dispersion 
of the f-fermions and inter- and intra- Coulomb repulsion. 



fermions compounds. 

We were successful enough to account for most of 
the experimental features. First, we have explained 
why there are seemingly two apparent QCPs which 
fit at the right respective coverage. The experimental 
one results from an extrapolation to zero temperature 
of an intermediate energy regime, while the theorit- 
ical one charachterizes the vanishing of the effective 
hybridization. We reproduced then the slopes and 
exponents of the coherence temperature and effective 
mass closely to the experimental results. The apparent 
lack of quantum critical behavior in temperature is 
qualitatively explained by the remarkably low energy 
scale of the spin liquid parameter on the ordered side, 
and the coincidence between the coherence temperature 
and the effective mass of the spinons in the hybridized 
one. Finally, we recover the fact that the activation 
gap, observed experimentally, has to vanish in the Fermi 
liquid phase before the critical coverage is reached, right 
when the system enters a strong hybridization regime 
for which the upper hybridized band becomes empty. 

Our study suffers though from some weakness and 
drawbacks. We used 4 fitting parameters, 3 for the hy- 
bridization and one for the slope of the coherence tem- 
perature. This is expected in any mean-field appraoch, in 
particular owing to the crucial role of the hybridization 
for the Kondo breakdown QCP, and can not be avoided 
at this level. The fact that the number of holons is 
too big away from the QCP, especially near the promo- 
tion coverage of the second layer, restricts the domain 
of validity of our model very close to the QCP. Finally, 
magnetism on the ordered side of the phase diagram is 
not handled in our model. Magnetism is best considered 
in the so-called slave fermions approach, which in turn 
describes badly the hybridized phase. 

But still the model is simple and strong enough to 
make predictions and put them to the test. This is the 
only proposed model of a de-confined QCP to be tested 
from ab-initio parameters. 

Useful discussions with H. Godfrin, G. Misguich, M. 
Neumann, J. Nyeki, O. Par collet, M. Ferrero and J. 
Saunders are aknowledged. 

This work is supported by the French National Grant 
ANR26ECCEZZZ. 



We benefited from the extensive literature on He^ to 
extract carefully most of the parameters of the model 
from the bare parameters. Crucial parameters, like the 
hybridization, were used as fitting parameters owing to 
the level of approximation of our study. Finally, we have 
emphasized some differences with intermetallic heavy 
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APPENDIX A: EVALUATION OF SOME 
INTEGRALS 

In here, we will evaluate the integrals in the mean- 
field equations (25). At T = 0, the calculation of these 



integrals is analytical for linearized bands in which case 



Po I de, 

D 



(Al) 



where po is the density of states at the fermi surface. 
Let's call 

T ^ G//(k,iu;„), 



A 
B 
C 

V 

£ 



T ^ Pj,(k,iw„), 



kcr,c 



T ^ ekG//(k,icjn) , 

k,cr,a;ri, 

dGff 



T ^ aek- 



dh 



db 



We diagonalize the 2x2 matrix which accounts for the 
hybridization of the f- and c- bands: 



-Ski = 



1 



A = (e°-ek)2+4(5y)2 . 

The integrals are all performed in the same way, first by 
summing over the Matsubara frequencies , and second 
by doing the momentum integration. The momentum 
integration is done by linearization of the band. 



A = 2TY^ 



kjUJn 

rD 



2p. 



dz 



2i7T {z-E_){z-E^) 
where the contour is on the whole complex plane 

nF{E_){E_ - Ck) nF{E^){E^ - e^) 

lit I 

D \ 



2po I de 



po 



po 



de- 



de- 



V (E- - E+) 
-y + ^y-' + AjbVY 



{E- - E+) 



-y - + A{bV f 
^y-^+A{bVY 



with Cm and the Fermi levels for the upper and lower 
bands respectively. 



(-6/ + - ^(e/ + a'^Y + Aa' (hVf) / {2a') 



ep = {-Cf + -^^{cf^ a'^Y + Aa' {hVY) / (2a') 

with the conditions —D < < 0; < < I) , 
and a' = alP' -\- <po/D . 

One obtains 

A = [-2y_D +ym- Vvl + 



ym. = {l- a')<^m -ef - i^; 
where j/p = (1 — a')€p — ef — fi; 

y-D = -(1 - a')D -ej - IJ . 

We proceed in the same way for B, C, V and 8 to find: 



B 



2po 



(1 - a') 



In 



Vm + Vyl + ^bvy 

yp + ^yl + A{bV f _ 



\yp + ^yl + i{bvr 
+ (ef + i^)ym + yl, 

- {y^/2 + ef + ^)^yl+4ibVr 

+ {^f + ij)yp + yl 

+ {yp/2 + ef + y)^yl + A{bVf 



V 



iaV^bpo 
(l-a')2 



In 



ym + Vyl + ^bvy 
yp + Jyl + 







(!-«') 









+ 



2V%^^y^ + W%^ 
-2 



(yn. - ./yl + W%^)^yl + W%-' 
-2 

{-yp - ^yl+W%^)^yl+AV%\ 
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APPENDIX B: CORRECTIONS TO SCALING 
FOR THE HOLON MASS 

In this appendix, we discuss the stability of the QCP. 



We will start by evaluating the diagram Ua 



f 



(1) 



f 



n,m^O k,q 

Introducing the angle 6 defined by ek+q = ek + v^qcosO 
and considering linearized bands like in (Al), we have, 
with g = SPpV^/po 



Ode- 



2 I l\^m\ 

^ a'q 



1 



1 



{icOn + i^m — — a^VFQCOsO — Cf) 

Summing over the fermionic Matsubara frequencies oOn 
then integrating over e, we get 



ni'^ = gpoa'T ^ / d0 



q,m7^0 
1 



Q 



2 ^ 7\^ru\ 

Oi' q 



1 



1 



i^m — a'vFQCosO — a' 11 — e/ 
Now, we have q = (qx^Qy) with = cos^ and qy 



qsinO. We suppose qy qx and expand q = ^ q^ + q^ 
We find 

J|g,|+gg/(2|g,|)) 

1 

X — 

ia'yljn — a'vFqx — ~ ^/ 
1 



where A is an ultra-violet cut-off. 

A logarthmic singularity in q^ arises when we integrate 
over qy^ and keeping only this singular part, we write 



nil) ^ ^ 



87r2 



ql Log(A/|g^|) 

{iQrn - C^'VFqx) 



3 • 



ni^^ is performed by continuation in the upper half plane, 
if < and in the lower half plane if 0.^ > so that 
to avoid the pole in the Green's function (See Fig. 14). 




FIG. 14: Gontour of integration : the cross stands for a pole 
and the hatched line for a branch cut. 

Changing variables in q^ = iz we get 

^(1) _SPO^rp io^'^m 



87r2 

Z7rSgn(r^^ 
{a'vF)^ 



H)3(|l], 



■ a'vFz) 



Log 



A 



To perform the summation over m, we notice that 
T Yl!^\/T 1 = 2A is independent of T. The same sum 
without the m = term will be 2A — T and to logarith- 
mic accuracy, we obtain : 



T^Log 



A 



-TLog 



A 



where the dots stand for 0{T) terms. 
Finally 



nW(r) 



(Bl) 
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This term is of negative sign and dominant compared 
to thus it can destabihze the regime. It therefore 
puts the intermediate regime in a fragile situation. This 
is due to the presence of the fermion loop gp- 



c 



94 



/t 9F \f 



and 




FIG. 15: The fermionic loop and the g4 vertex 



Fortunately, in D = 2, a mode-mode coupling con- 
stant ^4, coming for exemple from the term —JniUj/A in 
Eqn.([5|, provides corrections to scaling of the same tem- 
perature dependance but with a positive signe —TlogT^ 
competing with the one calculated previously. Indeed, 
for ^4 > 0, the (/)^-theory is stable and, close to a QCP, 
the corrections to scaling follow the law^ 



mb{T) 



b b 



rpid+z-2)/z 



Precisely, the leading logarithmic contribution coming 
from the 54 vertex reads 



mb(T) 



67r ^ V T 



2D ^ 



(B2) 



The stability of the intermediate regime is then a mat- 
ter of prefactors between the two terms. It can lie on a 
fragile basis, as it requires strong enough ferromagnetic 
short range fluctuations. However, it has been shown 
that this regime is stable for 1) = sP, we can thus ex- 
pect that a small three-dimensional character could cure 
this instability. The correction to the boson mass coming 
from the gauge fluctuation 



^A/^^yvvv^)v\/^ 



APPENDIX C: EXPRESSION FOR THE FREE 
ENERGY 



We start with the Luttinger-Ward formula (31). Fmf 



is the sum of the free energies of the two bands, given by: 

± n 

The sum over the fermionic Matsubara frequencies is for- 
mally divergent but its temperature dependence can be 
extracted using the following spectral representation 



1 



doc 



TT I X — lUOr) 



then, performing the summation over Matsubara fre- 
quencies, with 

T > : = - - nF{uJn)^ 



we get 



We turn now to the bosonic part of (31 ) 



Fh 



E 



T 
47r 



(2^)2 log [D-H^,^m)] 

r+oo 

^ / qdqlog 



a'q 



The integral over the holon momentum is dominated by 
large momenta, and we have 



+ CXO 



m 

8a' 



7^ 



Summation over the bosonic Matsubara frequencies is 
performed in the same way as for the sum over fermionic 
frequencies, and we find, for the T-dependent part of it 



3 



We end up with the total free energy given by 



F = -— Mp.^P-)^^ 



(CI) 



soes like in D = 2 and are sub-dominantP^. 



17 



M. Neumann et al, Science 317, 1356 (2007) 

G. Stewart, Rev. Mod. Phys. 56, 755 (1984); 73, 797 
(2001). 

H. V. Loneysen et al. Rev. Mod. Phys. 79, 1015 (2007). 
P. Coleman et a/., J. Phys. Cond. Matter 13 R723 (2001). 
A. Rosch et a/., Phys. Rev. Lett. 79, 159 (1997); A. Rosch 
ibid 82, 4280 (1999). 

M. Roger, Phys. Rev. Lett. 64, 297 (1990). 

G. Misguich, B. Bernu, C. LhuilUer and C. Waldtmann, 
Phys. Rev. Lett 81, 1098 (1998). 

H. Franco, R. E. Rapp and H. Godfrin, Phys. Rev. Lett. 
57, 1161 (1986). 

D. S. Greywall, Phys. Rev. B 41, 1842 (1990). 

E. Cohin et al Phys. Rev. Lett. 86, 2447 (2001). 
K.-D. Morhard, C. Bauerle, J. Bossy, Y. Bunkov, S.N. 
Fisher and H. Godfrin Phys. Rev. B 53 , 2658 (1996). 
M. Roger, C. Bauerle, H. Godfrin, L. Godfrin and J. 
Treiner, J. Low Temp. Phys. 112, 451 (1998). 
A. Hewson, in The Kondo problem to heavy fermions, 
Cambridge U. Press, page 333. 

A. Benlagra and C. Pepin, Phys. Rev. Lett. 100, 176401 
(2008). 

L Paul, C. Pepin and M.R. Norman, Phys. Rev. Lett. 98, 
026402 (2007). 

C. Pepin, Phys. Rev. Lett. 98, 206401 (2007). 

C. Pepin, Phys. Rev. B, 77, 245129 (2008). 

D. Vollhardt, Rev. Mod. Phys.56, 99 (1984). 

M. Heritier, Le journal de Physique - Lettres, 40, L-451 
(1979). 

^° P. Coleman, Phys. Rev. B 29, 3035 (1984). 
2^ S. Tasaki, Prog. Theor. Phys. 79, 1311 (1988). 

L. Pricaupenko and J. Treiner, Phys. Rev. Lett 72,2215 

(1994). 

S. Burdin, A. Georges, and D. R. Grempel, Phys. Rev. 
Lett. 85, 1048 (2000) 

P. Nozieres, Ann. Phys. Fr. 10, 19-35(1985); P. Nozieres, 
Eur. Phys. B 6, 447-457(1998). 

LPaul, C. Pepin and M.R. Norman, Phys Rev B, 78, 
035109 (2008). 

A dynamical exponent z is defined as the number of space 
dimensions that are effectively taken by the quantum fluc- 
tuations, namely u ^ . 

A.V. Chubukov , D. L. Maslov, S. Gangadharaiah, L. L 
Glazman Phys. Rev. B 71, 205112 (2005). 
The fitting of the data with the three parameters depends 
mainly on adjusting the hybridization and is very robust 
to a change in the initial condition, which can always be 
accounted by a rectification of the parameters Vi and V2. 
Taking into account the self- consistent determination of b 
changes very slightly the result of Refi^ and the phase di- 
agram is essentially the same. 

2^ P. A. Lee, N. Nagaosa and X-G. Wen, Rev. Mod. Phys. 
78, 17-85 (2006). 

^° D. Belitz, T.R. Kirkpatrick and T. Vojta, Phys. Rev. B 
55,9452 (1997). 

J. Zinn- Justin, Quantum field theory and critical phenom- 
ena, Oxford Science Publication. 

J. Rech, C. Pepin and A. V. Chubukov, Phys. Rev. B 74, 
195126 (2006). 



